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r— i Abstract 

-5: 

^ \ A complete calculation of the corrections to pion-deuteron scattering 

length up to order x 3//2 with \ — m -w/M N is performed. The calculation in- 
cludes the dispersive contributions and corrections due to the explicit treat- 
ment of the A resonance, s-wave pion- nucleon scattering parameters are 
extracted from a combined analysis of modern experimental data. 

cn : 

^ ! 1 Introduction 

The pion-nucleon (ttN) scattering lengths are fundamental quantities of low- 
energy hadron physics since they test the QCD symmetries and the pattern 
of chiral symmetry breaking. As stressed by Weinberg long time ago, chi- 
.£h ! ral symmetry suppresses the isoscalar nN scattering length a + substantially 

compared to its isovector counterpart a~ . Thus, a precise determination of 
a + demands in general high accuracy experiments. 

Here pion-deuteron (ttcL) scattering near threshold plays an exceptional 
role for Re(a 7r( i) = 2a + + (few-body corrections). The first term ~ a + is 
simply generated from the impulse approximation (scattering off the proton 
and off the neutron) and is independent of the deuteron structure. Thus, 
if one is able to calculate the few-body corrections in a controlled way, ird 
scattering is a prime reaction to extract a + (most effectively in combination 
with pionic hydrogen measurements). 
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A method how to calculate processes on few nucleon systems with external 
probes was proposed by Weinberg in one of his classical papers [1]. As a first 
step the perturbative transition operators need to be calculated using the 
rules of ChPT. Then those transition operators must be convoluted with the 
appropriate NN wave functions. This scheme was already applied to a large 
number of reactions like ird — > red [2], 7c? — > ir°d [3,4], 7r 3 He— > 7r 3 He [5], 
ix~d — > 7nn [6], and jd — > -n + nn [7], where only the most recent references 
are given. The standard expansion parameter x — m n/MN, where m n (Mjv) 
is the pion (nucleon) mass, was used in most of these references. 

It was also Weinberg who calculated the leading order few body correc- 
tions to the Trd system, the most important of which - the diagram when pion 
rescatters on two nucleons with the Weinberg- Tomozawa (WT) vertices - is 
almost as large as the experimental value for the ird scattering length. The 
diagrams at leading order were calculated in the fixed center kinematics, i.e. 
with static nucleons. An accounting for the nucleon recoils leads to poten- 
tially sizable corrections of order of % 1//2 in the standard Weinberg counting. 
The non-analyticity of this correction is related to the few-body singularities 
that are employed in some pion-few-nucleon diagrams as demonstrated in 
Refs. [7,8]. It was the main result of Ref. [8] that the importance of the 
resulting effect of all recoil terms is directly connected to the Pauli principle 
for the nucleons in the intermediate state. In particular, if the s-wave NN- 
state is not allowed by quantum numbers, which is fulfilled in the Trd process, 
the net effect of the recoil correction is to be small due to a cancellation of 
individually large terms. At next-to-leading (NLO) order there are basically 
the same diagrams as at LO but with subleading vertices. The calculation 
performed in Ref. [2] showed that the sum of diagrams contributing at NLO 
vanishes. Furthermore, the solution for {a + , a"} was found in Ref. [2] from a 
common intersection of three bands corresponding to the shift and width of 
pionic hydrogen atom [9] and to the shift of pionic deuterium [10]. Due to a 
cancellation of terms at orders x 1 / 2 and x the results for {a + ,a~} extracted 
in ChPT [2] turned out to be quite similar to the phenomenological calcu- 
lations [11,12]. However, the recent measurement of the width of n~p atom 
with much better accuracy [13] seriously changes the picture. The problem 
is that an intersection region of the three bands and thus a unique solution 
for {a + ,a~} docs not exist anymore. This finding means that something 
important is missing in our understanding of the nd system. This could be 
isospin symmetry breaking (ISB) effects that were recently found in Ref. [14] 
to give a huge effect to the nd scattering although with large uncertainty. 
Another possibility is that higher order effects to the transition operators 
could be important. In this presentation we discuss both possibilities. The 
influence of ISB effects on the extraction of the s-wave ttN scattering lengths 
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is considered in sec. 2. In sec. 3 and 4 we discuss corrections to the ird 
scattering length emerging at order N 3//2 LO (x 3 ^ 2 )- At this order two new 
classes of diagrams start to contribute. One of them, the so-called dispersive 
correction due to the process ird — > NN — > ird is the subject of sec. 3. In 
sec. 4 we discuss the effect of the A isobar as explicit degree of freedom. The 
main results are summarized in sec. 5. 

2 ISB effects and s-wave nN scattering lengths 

Since the leading one-body contribution (~ a + ) to the ird scattering length is 
chirally suppressed the role of ISB effects in this process becomes significant. 
For the ir~d system so far only leading ISB corrections were evaluated [14]. 
They were found to give a very large effect of order of 40% to the ir~d 
scattering length. In this section we would like to investigate the influence 
of this correction on the combined analysis of experimental data and thus 
on the s-wave ttN scattering lengths. To account for the ISB correction at 
leading order we should replace 2a + by a n - p + a n - n in the expression for a nd , 
which agrees to the former only, if isospin were an exact symmetry. The 
expressions for the ttN amplitudes with inclusion of ISB effects were derived 
in Ref. [14]: 



' 4(m 2 - m 2 (1 ) _ e 2 

X ) V n Cl ~ 2 
'4(m 2 -m 2 ) e 2 



a-Tr-p = a + a + A _ n ■ — [ — c x - — (4/i + h) 



Here m 7r (m 7r o) is the charged (neutral) pion mass, F w = 92.4 MeV and c\ 
and /i,/a correspond to the strong and electromagnetic LECs. Whereas 
f2 and ci are known more or less well (/2 = —(0.97 ± 0.38) GeV -1 [15] 
and ci = — 0.9±g;| GeV" 1 [16]) the value for h (\fi\ < IA GeV" 1 ) is very 
uncertain - a naive dimensional analysis was used in Ref. [14] to fix the 
latter. At the same time it is fx and Cx that give the largest contribution to 
the ISB correction for nd scattering thus introducing a large uncertainty in 
the extraction of {a + , a - } from the data. At this stage we would like to note 
that the parameters a + , cx and fx enter the expressions for a^- p and a n - n in 
the same linear combination (see Eq. (pQ)). Note that the expression for the 
charge exchange amplitude ir~p — > ir°n does not depend on the LECs c\ and 
fx at all. Therefore let us introduce the quantity d + which is defined as 

~+ + , 1 / ^(m 2 -m 2 ) 2 \ 

" + MTTx){ Fl <x-2efx). (2) 
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The leading isospin breaking terms that are also the main source of the 
uncertainty 2 ] are contained now in a + . Using this the expressions for a w - p , 
a^-n and Re a n d take the form 



a n - p — a + + a — , — r— f'2, 



1 e_ 2 

4tt(1 + X ) 2 



1 e 2 

O'n-n — a —a -+- - — — — ■ — -— j 2 , 

4vr(l + x) 2 

Kea n d = 2a + + (few-body corrections (a - )) , (3) 

Thus, we get a system of three equations for a n ~ p , Re a^d and a 7r - p ^ 7r o n (the 
explicit expression for the latter is given in Ref. [14]) to determine a + and 
a~ . Once this determination is performed and provided that new and less 
uncertain information about the LECs c\ and j\ is available from elsewhere, 
one will be able to extract a + directly without doing the analysis of pionic 
data once again. Let us discuss Eq. (j3J) in more detail. Note that the equa- 
tions for the hydrogen and deuterium shifts (a n - p and Rea^ respectively) 
written in terms of a + and a~ and including ISB effects at leading order are 
very close to those obtained in the isospin symmetric case in Ref. [2] in terms 
of a + and a". The difference basically consists in the term proportional to 
/ 2 for the pionic hydrogen shift that gives a relatively small effect. Thus 
the main modification due to the inclusion of ISB effects at leading order 
consists in the replacement of a + by a + . At the same time we know that 
for the isospin symmetric case there is no unique solution for {a + , a - } if the 
new data for the hydrogen width are utilized. Thus we conclude that the 
system of equations for {a + ,a~} does not have a unique solution either at 
least as long as the few body corrections of Ref. [2] are used. In the following 
sections we present the recent progress in this sector. 



3 Dispersive corrections 

Experimental measurement of the ird scattering length shows that its imag- 
inary part is relatively large, about 1/4 of its real part [10]. The imaginary 
part can be expressed in terms of the ird total cross section through the 
optical theorem. One gets 

47rlm(a 7rd ) = lim q {a(nd -> NN) + a(icd -> 7iViV)} , (4) 

2 The idea of using some linear combinations of observables to reduce the uncertainty 
was suggested in Refs [14, 17]. In particular the combination 2a 7r - p — a^- d that depends 
solely on a~ was considered in Ref. [17]. 
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a) b) c) d) 



Figure 1: Dispersive corrections to the ird scattering length. 

where q denotes the relative momentum of the initial ird pair. The ratio 
R = lim^o {a {ird -> NN)/o(ixd -> jNN)) was measured to be 2.83 ± 0.04 
[18]. At low energies diagrams that lead to a sizable imaginary part of some 
amplitude are expected to contribute also significantly to its real part. Those 
contributions are called dispersive corrections. As a first estimate Bruckner 
speculated that the real and imaginary part of these contributions should 
be of the same order of magnitude [19]. This expectation was confirmed 
within Faddeev calculations in Refs. [20]. Here we present the first consistent 
ChPT calculation of the dispersive corrections that emerge from the hadronic 
ird — > NN — > ird and photonic ird — > 'yNN — > ird processes [21]. We define 
dispersive corrections as contributions from diagrams with an intermediate 
state that contains only nucleons, photons and at most real pions. Therefore, 
all the diagrams shown in Fig. [1] are included in our work. All these diagrams 
contribute at order x 3//2 as compared to the leading double scattering dia- 
gram(see Ref. [21] for details). The hatched blocks in the diagrams of Fig. [1] 
refer to the relevant transition operators for the reaction NN — > NNtt that 
consist of the direct and rescattering mechanisms. Note that the latter is to 
be calculated with the on-shell ttN — > ttN vertices (2m n ) as was derived in 
Ref. [22]. Using the CCF potential [23] for the NN distortions we found for 
the dispersive correction from the purely hadronic transition 

Sa^ p = (-6.5 + 1.3 + 2.4 - 0.2) x 10~ 3 m; 1 = -3.0 x 10" 3 m~ l , (5) 

where the numbers in the first bracket are the individual results for the 
diagrams shown in Fig. [TJ in order. Note that the diagrams with interme- 
diate NN interactions and the crossed ones (diagram (c) and (d)), neither 
of them were included in most of the previous calculations, give significant 
contributions. When repeating the calculation with the four different phe- 
nomenological iVW potentials CD Bonn [24], Paris [25], AV18 [26] we find 

5a*f = (-2.9 ± 1.4) x 10- 3 m- 1 , (6) 

where the first number is the mean value for the various potentials and the 
second number reflects the theoretical uncertainty of this calculation esti- 
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mated conservatively — see Ref. [21] for details. Note that the same calcu- 
lation gave very nice agreement for the corresponding imaginary part [21]. 

In Ref. [21] also the electromagnetic contribution to the dispersive cor- 
rection was calculated. It turned out that the contribution to the real part 

was tiny 0.1 x lO^m^ 1 — while the sizable experimental value for the 

imaginary part was described well. 

4 Role of the Delta resonance 

From phenomenological studies it is well known that the delta isobar A (1232) 
plays a very special role in low energy nuclear dynamics [27] as a consequence 
of the relatively large 7riVA coupling and the quite small delta-nucleon mass 
difference A = Ma — Mn ~ 2m n , where Ma denotes the mass of the delta. 
In the present section we investigate the role of the A isobar in the reaction 
red — > red at threshold in EFT. In the delta-less theory the effect of the A 
resonance is hidden in the LEC c<i which is the leading term in the so-called 
boost correction to the ltd scattering length [2]. This correction is known 
to be quite sizable (~ 10-20% of a n d) although very model dependent. The 
pertinent one-body operator scales with the square of the nucleon momen- 
tum and therefore the corresponding expectation value is proportional to 
the nucleon kinetic energy inside the deuteron — this quantity is strongly 
model-dependent [28]. However, the value of C2 is reduced by a large factor 
once the delta contribution is taken out [29, 30] so that the residual boost 
correction becomes negligible [31]. 

The reason why the explicit inclusion of the delta in pionic reactions on 
the two-nucleon system is beneficial is that the dynamical treatment of the 
A allows to improve the convergence of the transition operators. Let us, 
for example, focus on the one-body terms with the delta (see, e.g., second 
diagram in Fig. EJ). Then the corresponding tcN — > ttN transition potential 
is proportional to p 2 / (m n — A — p 2 /M^). For static deltas, the nucleon-delta 
propagator reduces to l/(m n — A). Thus, in the latter case the transition 
operator behaves like p 2 , whereas in the former it approaches a constant for 
momenta larger than \p& | ~ a/(A — m^M^ ~ 2.7m n with the effect that 
the static amplitude is much more sensitive to the short range part of the 
deuteron wave function and must be balanced by appropriate counter terms. 
The value of pa is numerically very close to p t hr = ^/M^ra^ ~ 2.6m w - 
the minimum initial momentum for the reaction NN — > NNtt (for a recent 
review of this class of reactions see Ref. [32]). This automatically puts the 
delta contributions in the same order as the dispersive corrections [31]. In 
Fig. [2] we show diagrams with the A isobar that contribute at order x 3 ^ 2 - 
Diagrams with crossed external pions are not shown explicitly but are taken 
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Figure 2: Diagrams with the A resonance at order x 3 ^ 2 ■ Crossed terms are 
not shown explicitly but included in the calculation (see Ref. [31]). 



into account in the calculation. The resulting correction is [31] 

5a* d = (2.38 ± 0.40) x 10" 3 m' 1 , (7) 

where the central value is the arithmetic average of the results for the seven 
different potentials and the uncertainty reflects the variations in the results. 
Note that we used phenomenological NN models without [24-26] and with 
[23] explicit delta degree of freedom, as well as three variants of NN wave 
functions derived within EFT [33]. These numbers were obtained with the 
7r N A coupling constant Ha — 2.77. In contrast to earlier treatments of the 
boost correction, the results we found with the explicit treatment of the A 
depend only very weakly on the NN model used. In Ref. [31] also a detailed 
comparison to previous phenomenological works is given. 

5 Results and Conclusions 

We performed a complete calculation of the isospin-conserving corrections 
to the pion-deuteron scattering length up to order x 3 ^ 2 - The calculation 
includes the dispersive contributions and corrections due to the dynamical 
treatment of the A resonance. Although these corrections are quite signifi- 
cant individually the net effect of the diagrams that contribute at order x 3 ^ 2 
is very small: 

+ = ("0-6 ± 1.5) x 10- 3 m-\ (8) 

However, an important consequence of our investigations is that once the 
A is treated dynamically, as it is done here, the so-called boost corrections 
contribute insignificantly to the ird scattering length. 

Also we analyzed the role of ISB effects at leading order in the combined 
analysis of pionic data. It was observed that the LEC fx, that is known 
very poorly, appears in the expressions for a v - p and a^- n in the same linear 
combination with a + and the LEC C\. We called it d + . Thus, the inclusion 
of ISB effects at leading order consists basically in the replacement of a + 
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by a + in the combined analysis of pionic data. This drastically reduces 
the uncertainty of the analysis that originates mainly from our ignorance 
regarding f\. The solution for the s-wave ttN parameters {a + ,a~} is shown 
in Fig. [3j The black band stems from the analysis of the pionic hydrogen shift 
[13]. The blue vertical band corresponds to the new preliminary data for the 
pionic hydrogen width [34]. The red solid and dashed bands correspond to the 
pion deuteron scattering length calculated with and without corrections at 
order x 3//2 - Also the boost correction was not included in the full calculation 
corresponding to the solid red band as a consequence of the explicit treatment 
of the A resonance. It is basically the latter effect that improves the situation 
resulting in some intersection region for all three bands. However, it still 
remains to be seen if the corrections at NLO of the isospin violation do 
not distort this picture. Corrections at this order for the 7r~p system were 
evaluated in Refs. [15,35] and turned out to be quite sizable, especially those 
that come from the pion mass difference. In order to push also the calculation 
for the ird system to a similar level of accuracy in isospin violation, the 7r~n 
scattering amplitude as well as some virtual photon exchanges in the n~d 
system are still to be calculated. 
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